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Abstract. Differential equations are derived for a continous limit 
, of iterated Schwarzian reflection of analytic curves, and solutions 

\& ' are interpreted as geodesies in an infinite-dimensional symmetric 

space geometry. 

O: 

Q; 

ri ■ 1. Introduction 

^ . 

Given a pair of (sufficiently nearby) regular analytic curves Ti and 
T 2 in the complex plane, Schwarzian reflection of 1^ in T 2 produces a 
third such curve r 3 = T 2 • r x . Iterating this construction yields curve- 
dynamics {T n+1 = r ra -r n _!} with discrete-time n. Letting T 2 approach 
qq \ T\ and taking a second order approximation of the above construction 

(N ■ yields a quadratic second order PDE (Equation [32|), which integrates 

^ . at once to a first order linear PDE (Equation p0|) . 

O '. The resulting continuous reflection process {T t } describing confor- 

mally invariant evolution on A = {T} (unparametrized analytic curves) 
has a compelling interpretation in terms of infinite-dimensional geom- 
etry on A. We observe that Schwarzian reflection of curves satisfies the 
formal properties of a symmetric space multiplication //(r 2 , Ti) = r 2 Ti 
(Theorem [I]), as defined in ||Loos]| . Further, the canonical affine con- 
nection V on A determined by fi yields precisely Equation |32| as the 
^ . geodesic equation on A (Theorem |3]). 

Our formulation makes essential use of the Schwarz function S(z) 
as natural parameter on A. Following |Davisj ], the Schwarz function of 



O 

S 



a simple analytic arc T is a holomorphic function defined near T such 
that T is implicitly defined by T = {z : z = S(z)}. (For non-simple T, 
S(z) is defined on a Riemann surface.) The subsequent identification 
of points in A with Schwarz functions A = {T} ~ {S(z)} yields a calcu- 
lus perfectly suited to the symmetric space structure on A — in fact the 
Schwarz functions play the role of transvections — and our PDEs govern 
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time-evolution of Schwarz functions S(t, z) representing geodesies. Al- 
ternatively, the first order PDE describes one-parameter subgroups in 
the image of the Cartan immersion C : G/H — > G, and left /i-translates 
of same. (In the standard- looking identification A ~ G/H, we note G 
and G/H have only formal meaning, explained below.) 

The resulting geodesic equation is trivial to solve. When the first 
order PDE is reduced to an ODE (say, by the method of characteris- 
tics), one finds that the normal velocity field along the moving curve 
Ft is just (the real part of) the holomorphic extension W(z) = iv(z)dz 
of the velocity field along the initial curve IV As an orthogonal pair of 
(singular) foliations, the curves F t and the trajectories of their moving 
points are nothing but the equipotentials Re[u] = and streamlines 
Im[u] = of a meromorphic differential u = dz/iv(z) (Theorem [|). 

We contrast the motion of points on F t with the classical dynam- 
ics of particles in stationary ideal flows (see, e.g., ||Cohn|| , |Springcrj| ) . 



Whereas W is dual to u — as vector to covector — the duality associat- 
ing uj and a fluid velocity field V (the gradient of a harmonic function) 
depends in fact on a metric ( , ), and is reciprocal to W; that is, 

V = W/(W,W) (for any conformal ( , )). Thus, e.g., the replacement 

V — > W, turns a vortex point for the fluid into a stationary point about 
which F t pivots, a source into an expanding closed curve, and a dipole 
into a stationary point of first order contact for the family F t . 

In particular, these three behaviors are exhibited by spacelike, time- 
like, and lightlike geodesies (respectively) in the three-dimensional sym- 
metric subspace A 3 C A consisting of circles in the Riemann sphere, 
with multiplication \x given by circle inversion (Examples [T^). As 
the above nomenclature suggests, A 3 C A turns out to be Lorentzian, 
with Levi-Civita connection coinciding with the canonical symmetric 
space connection V. One thus recovers Mobius circle geometry (see 
|(Jecil|| ) , and the geodesies S(t,z) are the elliptic, hyperbolic, and para- 
bolic pencils of circles of classical geometry. 

Away from stationary points, all geodesies are locally conformally 
equivalent; at a stationary point zo, the tangent line to F t has a well- 
defined constant rotation rate 9{zq); when 8 = 0, the Kasner horn angle 
invariant JC = 2 /O s changes at a constant rate (Proposition [TJ). The 
order of a pole and the residue of uj = dz/iv(z) distinguish locally 
inequivalent geodesies, while periods of uj carry global information. 

Other significant features arise if singularities are allowed in the an- 
alytic curves themselves. For instance, confocal ellipses will be seen 
to describe geodesies (Example [14]), in which continuation of S(t,z) 
through the time t at which branch points are encountered (at the 
common focii) pushes the reflection process through an apparently 
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singular curve T to (the branch cut). The local behavior at branch 
points — where T t is stationary only instantaneously — is easily under- 
stood in terms of an alternative representation of geodesies given by a 
quartic second order PDE (Equation [37]) for an evolving parametrized 
curve w(t, x). The latter equation admits special solutions of the form 
w(t,x) = f(x + ict), with f(z) holomorphic. In particular, rational 
functions f(z) yield a variety of simple examples of geodesies with sin- 
gularities given by f'(z) = — see Figures |3], f|. 

The equations of continuous reflection, the various phenomena exhib- 
ited by concrete examples, and the formal symmetric space structure of 
A are developed here as elements of Schwarz reflection geometry. In the 
same spirit, the sequel to the present paper ( ||Carini-Langer|| ) begins to 
address, among other topics: metric structures on A and relation to the 
Witt algebra; invariants/normal forms for tangent vectors under con- 
formal equivalence via the isotropy representation of G/H; holonomy 
of V, the group of displacements of A and the image of the exponen- 
tial map; finite-dimensional geometries embedded in A (defined, e.g., 
by constraints on singularities of S(z)); relation to theory of quadratic 
differentials and meromorphic flows ( |Strebe]f| , ||Mucino-Raymundo|| ) . 

To place our approach in context, we note that the general theory 
of symmetric spaces (like that of Lie groups) has only limited appli- 
cability to infinite-dimensional examples, which require individual at- 
tention. The manifold with multiplication approach of Loos may offer 
advantages in some such cases, owing to its directness. In particular, 
it lends itself to a simple computational formalism which has proved 
very convenient for the study of A. On the other hand, for topological 
reasons, one does not know precisely how A ought to be regarded as an 
example of a (locally) symmetric space, nor what would be gained as 
a result. In our view, the question itself is premature until the inter- 
esting particulars of this example are more fully explored — hence the 
stepping-stone of Schwarz reflection geometry. 

Some of the issues may be more familiar in the homogeneous space 
context. Geometric and Lie algebraic considerations tempt one to iden- 
tify A (unparametrized curves) with the quotient G/H (parametrized 
curves modulo reparametrizations), where H = Diff(S l ) (analytic 
diffeomorphisms of the circle), and G = He ("analytic embeddings 
of S 1 "). Unfortunately, H has no complexification H c (however, see 
Neretin|| , [ |Segal|| ). This point is discussed, e.g., in [Pressley-Segal 



where the "near-compactness" of loop groups and their homogeneous 
spaces is highlighted as exceptional behavior in the infinite-dimensional 
context. One should also contrast A with the much-studied coad- 
joint orbits of the circle diffeomorphism group (or the Virasoro group). 
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For instance, standard infinite-dimensional analysis provides a mani- 
fold structure on M = Dif f{S l ) / S l , hence, a satisfactory topolog- 
ical setting for the well-known Kahler and symplectic structures on 
M (and for related spaces of projective and conformal structures). 
For a sampling of interesting recent work involving such structures 
and related equations, see ||Guieu-(Jvsienko] | , | JKirillovj | , |Michor-Katiu |. 



| W iegmann- Z a br o din] | , and references therein. Some of this work is ul- 
timately connected to our subject, as will be explained further in the 
sequel. 

The organization of the paper is as follows. In §2, we discuss reflec- 
tion of analytic curves in terms of Schwarz functions and symmetric 
space multiplication. In §3, we consider Davis iteration, its continu- 
ous limit, and the implied partial and ordinary differential equations 
for continuous reflection. In §4, the second order PDE for continuous 
reflection is interpreted as the geodesic equation on A. Though sym- 
metric space formalism underlies much of §2-§4 and appears explicitly 
in a few places, we have emphasized concrete constructions and direct 
derivations of functional and differential equations. A reader unfamil- 
iar with symmetric space theory may refer to Appendix I as needed for 
a brief review of the most relevant aspects of symmetric space formal- 
ism. (We note that our emphasis is non-standard; in particular, we have 
minimized the homogeneous space viewpoint). Appendix //describes 
the canonical connection on an abstract symmetric space (following 
[Loos ]), and completes the proof of Theorem 



Finally, it is a pleasure to thank Jerome Benveniste and David Singer 
for many stimulating and helpful conversations. 



2. The Schwarz function and reflection 

We begin by recalling how an analytic curve V determines a reflec- 
tion Rr(z), as well as a Schwarz function S(z) = Sr{z). For a more a 
complete discussion, we refer the reader to the excellent MAA mono- 
graph The Schwarz function and its applications by Phillip J. Davis 
Pavisjl , and also to the more recent work [|5hapiro| |, which emphasizes 



certain analytical issues. We have made minor adaptations in the for- 
mer exposition, to suit our geometric (symmetric space) point of view. 

Reflection in the real axis 1 C C is given by complex conjugation, 
z i — ^ z — R(z) = Rr(z). A symmetric subset of C is one which is 
preserved by R. The conjugate of a holomorphic function is defined by 
the formula / i— > cr(/) = ctr(/) = R^ 1 o f o R; here, if / is defined on a 
domain U, then a(f) is defined on the reflected domain U. A function 
is symmetric if it is preserved by a. Note that symmetric functions are 
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precisely those defined on symmetric domains preserving symmetry of 
all subsets. 

In the present context, we will generally follow Davis' very convenient 
(but slightly dangerous) bar notation for conjugation of an analytic 
function f(z) = a n (z - z ) n : 

CO 

(1) f(z) = a(f)(z) = W) = E ~ 



In the abstract setting, a will denote an involutive automorphism of a 
group, realized here in especially transparent notation: 

(2) / = /, 7^9 = fog 

Reflection in R has induced conjugation of analytic functions, which 
in turn will be seen to yield Schwarzian reflection in analytic curves. In 
general, an (unparametrized) analytic curve will be regarded as a subset 
of the complex plane C which is the image T = 7(T ) of a standard 
base curve Y under a non-singular holomorphic function 7 : U — > C 
defined in a neighborhood U of r o . Here it will usually suffice to take 
the real axis r o = IR as base, though for some purposes the unit circle 
S* 1 = {z G C : zz = 1} will be the natural choice. A given curve may 
also be regarded as an equivalence class T = [7] = / ~fH consisting of 
all reparametrizations of a given parametrized curve 7 : r o — > C via 
analytic diffeomorphisms h G H of T . 

We may sometimes be required to regard Y as belonging instead to a 
Riemann surface S, e.g., when Y is not simple. £ may be constructed 
by analytic continuation, so that 7 is not merely non-singular (7^ 7^ 0), 
but a diffeomorphism 7 : U — > W C S. With this understanding, one 
easily verifies the following 

Proposition 1 (Reflection). Given Y an analytic curve, there is a do- 
main W D T and a unique antiholomorphic extension to W , Rv{z), 
of the identity map on Y. In fact, Schwarzian reflection in Y is 
well-defined near Y by: 

(3) R T = 1 oRo 1 ~ 1 , r = 7(R) 

Schwarzian reflection is involutive and conformally invariant: = id, 
and for ip a conformal map defined near Y, R v (r) = (p o R r o ip^ 1 . 

Observe that Rr may now be used in place of R to define T-conjugation 
of functions, err, and T-symmetry of sets and functions — with identical 
properties — effectively replacing r o = K with new base curve Y. We 
note also the relation to (a simplified version of) the Schwarz symmetry 
principle: if f(z) is defined and regular on a domain D containing Y, 
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then f(z) may be extended holomorphically to the symmetrized do- 
main D U Rr(D) by the formula Rfrv) ° / ° Rr- 111 the special case 
/(r) = T, the holomorphic extension of / to Rr{D) is T-symmetric. 
Parallel to the above proposition, we have 

Proposition 2 (Schwarz functions). Given T an analytic curve, there 
is a domain W D T and a unique holomorphic extension to W , S(z) = 
S-p(z), of the map z i— > z onT. In fact, the Schwarz function ofY is 
well-defined near T by the formula: 

(4) S v = RoR r =707"!, r = 7 (R) 

Not only does T determine S-p, but S determines T by the equation: 

(5) T = {zeC: S{z) = z} 

The last equation suggests turning the relationship between T and 
S around, and writing T = T s ; e.g., the identity function S(z) = z im- 
plicitly defines the real axis, R = Ts, and S(z) — 1/z defines the circle 
S 1 = Ts- However, for a typical holomorphic function S(z) = f(z), 
Equation [| yields not a curve, but a discrete set of points. The one- 
to-one correspondence between unparametrized analytic curves and 
Schwarz functions is regarded here as fundamental. In the following 
proposition, we collect some formulas which will be useful in relating 
Schwarz functions S to the geometry of the corresponding curves Tg. 

Proposition 3 (S(z) and curve geometry). LetT = 7(R) have Schwarz 
function S(z). Then the derivative of S, restricted to T, is a function 
of unit modulus known as the "clinant": 

(6) S z : T - RP\ S z = ~ x / lx = e~ 2ie 

The clinant may be used to describe the following geometric quantities: 

i) The unit tangent, up to sign: 

(7) i = e^ = \j^fS z 

ii) The signed curvature ofT: 

(8) k = \s zz /{S z f' 2 

Hi) The derivative of curvature with respect to arclength along T: 

(9) k 9 = ^{S,z} = ^-[(S zz /S z ) z - \{S ZZ /S z f] 

Here, {S, z} is the Schwarzian derivative. 

iv) Kasner's conformal invariant for a (normalized) horn angle: 

(10) K = k 2 / k s = ^S zz / S Z {S, z} 
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Proof. For Equation [| write 7^ = ^e td , and differentiate 7 = £(7) 
with respect to a; to obtain = 7^/7^ = \ie~ %e j \ie %e = e~ 2td (so 
w = /xe* 61 and w are homogeneous coordinates for S z G MP 1 ). We 
note that the clinant may be thought of as the "complex slope" of 
T, since the equation for the tangent line to T at z may be written 
(z - Zq) = S z (z Q )(z - Zq). 

Next, taking x to be an arclength parameter along 7, differentiation 
of 5* 2 = e~ 2td with respect to x gives S zz e l6 = — 2ie~ 2 * £l ^|, and solving 
for k — 4| gives Equation |8]. Equation || follows similarly by straight- 
forward computation, giving Equation [H] as well. For simplicity, we 
have specialized the formula /C = O 2 /0 S , O = «i — k 2 by assuming the 
second of the two curves in first order contact has been conformally 
mapped to a straight line (see ||Kasnerj|, footnote 13 on the 'natural 



measure of the horn angle' M\2 = K, between two curves in first order 
contact). Conformal invariance of K, may be derived by (double) appli- 
cation of the "chain rule" for Schwarzian derivative of f(w) = f(g(z)): 

{f(g(z)),z} = (g z ) 2 {f(w),w} + {g(z),z} 

We recall that {f(w),w} = iff f(w) is a Mobius transformation, so 
this well-known formula generalizes Mobius invariance of { , }. □ 

From the abstract symmetric space point of view, the description of 
Schwarz functions as transvections S = 7 o 7- 1 (relative to the involu- 
tion ctk) has many formal consequences. We begin with some definitions 
and identities for some basic operations on Schwarz functions: 

Proposition 4 (Inversion, multiplication, and Hermitian conjugation). 

i) The local inverse of a Schwarz function S = Sr near T satisfies 

(11) s- 1 = s, 

ii) The product of Schwarz functions S = Sa and T = Sb, 

(12) ^■T = Sor 1 o5 = 5oToS, 

is the Schwarz function of the product of curves A-B = Ra{B) — the 
reflection of B in A. 

Hi) If ip is a conformal map defined near V = T$, then the Schwarz 
function of the image curve ip(T) is given by Hermitian conjugation: 

(13) \(<p){S) = <p ° S o p- 1 

iv) Multiplication of Schwarz functions is conformally invariant: 

(14) X( V )(S-T) = X( V )(S)-X^)(T) 
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Proof. All claims may be regarded as symmetric space formalism — see 
the appendix. We verify ii) here, as multiplication of Schwarz functions 
is fundamental. Thus, let A = a(R) and B = /3(R). Then 

S ■ T = a o a -1 o p o P' 1 o a o a -1 = (a o a -1 o p) o (a o a -1 o 

so 5-T is the Schwarz function of the curve aoa _1 o/3(R) — in particular, 
Schwarz functions are closed under the above multiplication. Now, 
R A = a o R o a' 1 = a o a' 1 o i?, so = a o a" 1 o R o /3(R) = 

a o or x o /3(R), the curve whose Schwarz function is S ■ T. □ 

Conformal invariance of the product operation is closely related to 
the third property ( "anticonformal invariance" ) in the following 

Theorem 1. Multiplication of Schwarz functions satisfies the for- 
mal properties of multiplication on an abstract symmetric space. Namely, 
letting P, Q and R denote Schwarz functions, the following hold when- 
ever the products are defined: 

1. p.p = p 

2. P- (P-Q) = Q 

3. P-{Q-R) = {P-Q)-{P-R) 

4. P-Q = Q implies P = Q, or P andQ are "orthogonal" : P'(zo) = 
—Q'(zo) at intersection points (zq = P(z ) = Q(zq)). 

Proof. Properties (l)-(3) are formal consequences of the representation 
(Equation |j) of Schwarz functions as transvections. To prove (4), sup- 
pose P ■ Q = Q and P(z ) = Q(z ). Then p- 1 o Q = Q- 1 o P, so 
the function F = Q' 1 o P is involutive with fixed point z . Differ- 
entiation of F(F(z)) = z at zq implies F'(z ) = ±1. If F'(z ) = —1 
then P'(zo) = —Q'(zo). On the other hand, if F'(z ) = 1 note that 
H(z) = z + F(z) defines a local diffeomorphism (H'(z ) = 2) satisfying 
H(F(z)) = H(z). Thus, F(z) = z near z , hence, P = Q. □ 

Remark 5. Axioms (l)-(4) of an abstract symmetric space — see Def- 
inition |22| — may be paraphrased: for each P, the left-multiplication 
operator P- is an involutive symmetric space automorphism with iso- 
lated fixed point P. In view of Equation Property (4) above says 
that if Q 7^ P is a fixed point of P-, then the curves Tq and Tp can 
only meet orthogonally. In a C l or stronger topology on the space of 
Schwarz functions, P ■ Q = Q cannot hold for Q ^ P near P. For 
either Tq intersects Tp orthogonally, or Tq does not intersect Tp, but 
contains reflected pairs of points on either side of Tp — this would be 
impossible for any connected Tq lying in a thin neighborhood of Tp. 
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Example 6. Unoriented circles in the Riemann sphere, A 3 C A 

az+b 
cz+d 



Among linear fractional transformations, M(z) = ffr^, 5 = ad—bc^ 



0, the Schwarz functions are of the form: 

{ad — bc)z+ {ah — ab) ujz + iB 



(15) S{z) = M o M~ (z) 



{cd — cd)z + {ad — be) iAz + uj ' 



where u G C and A, B e R, with det(S) = ujuj + AB = 55 > 0. 
The corresponding curves Ts = {z : z = S{z)} satisfy the quadratic 
equation Azz + i{ujz — ujz) — B = 0, which describes all circles in the 
Riemann sphere S 2 . Given two such Schwarz functions Sj = ^ z+lB3 



jZ-\-LOj 



the symmetric space product S 2 ■ Si agrees with the usual inversion of 
circle l\ in circle T 2 (as is readily verified using conformal invariance). 
Thus, the three-dimensional space of unoriented circles in S 2 is indeed 
a symmetric subspace, A 3 C A. 

Further, setting u = x 2 + ixi, A = x 4 — rr 3 , B = —{x^ + X4), we 
introduce the homogeneous coordinate x = {x±, x 2 , x 3 , x 4 ) on the space 
of circles. Scaling x by real A ^ does not change S{z) or T; the 
normalization det{S) = (x,x) = x\ + x\ + x\ — x\ = 1 allows us to 
identify A 3 with antipodal pairs p = ±x of points in the one-sheeted 
hyperboloid Qf{l) = {x e M 3 : (x,x) = 1} in the Lorentz space 
Rf. Multiplication by p is an isometry with respect to the induced 
pseudo-Riemannian metric on A 3 (likewise for multiplication by x in 
the symmetric space of oriented circles Qf(l) — see the appendix). 

Now let inverse stereographic projection n^ 1 {z) = {f^=, in^z) ' f+il) 
define the coordinate z on S 2 = {{xi, x 2 , x$) : x\ + x\ + x\ = 1}. Then 
the circle T a C S 2 with homogeneous coordinate a = {a±, a 2 , 03, 04) is: 

(16) r a = S 2 fl {aiXi + a 2 x 2 + a 3 x 3 = a 4 } 

This may be verified directly using the above quadratic equation for 
T a = IV In the context of Mobius geometry, S 2 is embedded in R 3 
as the intersection of the light cone (x, x) = with the hyperplane 
X4 — 1. If P^- is the hyperplane through the origin R 3 -orthogonal to 
a, then T a = S 2 (~) P^. (Alternatively, S 2 may be identified with the 
projectivized light cone and Y a with a projective subspace V^.) 

3. Davis iteration and its continuous limit 

Given nearby initial curves, T , r 1; a sequence of Schwarz functions 
may be defined, inductively, by the Davis iteration scheme for {t„}: 

(17) S n+2 = S n+ i ■ S n = S n+ i o S n o S n+ i 
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We note that the iterates S n are nothing but the symmetric space 
powers (Si) n of Si relative to the base point So. For simplicity, we will 
discuss formal properties of the iteration for < n < oo. 

To formulate the continuous analogue of the iteration scheme, it will 
be convenient to appeal to the following 



Lemma 7. Davis iteration, Equation \F§ is equivalent to the following 
identity with integer indices < j, k: 

(18) S j+k = Sj o S o S k , 

Proof. First suppose Equation |18] holds, with initial Schwarz functions 
So, Si. We make the inductive assumption that S k is a Schwarz func- 
tion whenever < k < n + 1. Then S n+ \ = S n o So o S± implies 
b n o 6 n+ i = b o bi and hence b n+2 = b n+1 ob Q ob 1 = b n+1 ob n o b n+1 . 
Thus, S n+ 2 is also a Schwarz function, and Equation |T7] holds. 

Conversely, assume Equation Suppose, for some M > 0, Equa- 
tion [18] holds whenever < j + k < M. (Note that the identity is 
obvious when either j or k is zero — in particular, when M = 0). Let 
1 < j' — j + 1, k' — k + 1 be integers satisfying j' + kl = M + 1. Then 

So o o ooooooooo_ 

j'+k' — bj + k+l° bj + k° bj + k+l — bjl O £>0O bk° bQO SjO bjO £>0O £>k> — 



Sj/ o Sq o S)-', so Equation [L8| for < j, k follows by induction. □ 



Now we simply replace the integer variables j, k in Equation [18| with 
continuous "time variables" , t,u G M — thus interpolating Davis iter- 
ation (a heuristic argument justifying this step is given in the next 
section). The resulting equation for a time-dependent Schwarz func- 
tion S t (z) = S(t,z) defines continuous Schwarzian reflection: 

(19) S t+U = S t o S Q o S u 

Next, denoting t-derivatives by dots, = J^, and ^-derivatives by 
primes, '=■§;, Equation |1^ yields 

S t +u = ^^*+« = Q^St ° S ° S u = (S t o S o S U )(S' o S U )S U . 
Evaluating at u — 0, and using S' Q o So = I/Sq, we obtain: 

(20) S(t,z) = ^^S'(t,z)=g(z)S'(t,z) 

In this PDE, the coefficient g(z) = S(0, z)/S'(0, z) is regarded as a 
given holomorphic function, to be further characterized below. 
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Proposition 8 (Conformal invariance). Let S t (z) satisfy Equation [7^, 
hence also Equation Under conformal mapping by <f(z), the trans- 
formed Schwarz function, S t = X((p)(St) = (p o S t o tp" 1 , again sat- 
isfies Equation [H], and Equation \2Q as well, with g(z) replaced by 
g{z) = ^^-\z))g^~\z)). 

In view of the proposition, it appears reasonable to consider nor- 
malizations of the above equations in which the initial curve has been 
conformally mapped to the real axis: 5(0, z) = z is the "initial po- 
sition" and g(z) = S(0,z) the "infinitesimal variation". In this case, 
(|19"1) reduces to the functional equation for iteration groups: 

(21) S t+U = S t o S u , S = Id 

Further, z = S t o S t (z) implies = S (z) + S (z), so g(z) = S (z) is 
imaginary on the real axis. Accordingly, we let v (z) be the continuation 
of a prescribed real analytic function v(x) = ^S'o(x), and consider the 
normalized initial value problem: 

(22) S(t, z) = -2iv(z)S'(t, z), 5(0, z) = z, 

Remarkably, the Schwarz function S(t, z) now plays a second role: 

Proposition 9 (Canonical parametrization). IfS(t,z) satisfies Equa- 
tion then T t = {z : z = S(t, z)} has parametrization 

(23) x ^ 7(t, x) = S(-t/2, x), x e R 

With real analytic function v(x) as above, 7(t, x) satisfies the PDE 

(24) j(t, x) — iv{x)^'(t, x), r y(0,x) — x, 

hence, j(t,x) describes normal motion ofT t : ( A f(t,x), r y'(t,x)) = 0. 

Proof. By Equation |2l|, we have S^ 1 = S t = S-f Thus, S(t,j(t, x)) = 
S(t, S(-t/2, x)) = S(t/2, x) = S(-t/2, x) = 7(t, x) = 7(t,x), so >y(t, x) 
lies on T t . (In terms of symmetric space formalism, the appearance 
of the factor —1/2 may be traced to the special form of the Cartan 
immersion on cosets of transvections.) 

Since ^(t,x) = §f (— t/2,x) ^ 0, a regular parametrization of T t is 



thus obtained, and Equation ^4] follows from Equation |22|. The last 
statement follows at once, and gives direct geometric meaning to the 
function v(x) as the initial rate of displacement of T t from r . (Similar 
interpretation of the extension v(z) is obtained below.) □ 
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By invoking the method of characteristics and subsequently applying 
the group property for S(t, z), one may reduce Equation [22| to the ODE: 

(25) ^ = -2iv{S), S{0,z) = z 

A more direct derivation is contained in the proof of the following 

Theorem 2. The equations of continuous reflection are reducible to 
quadrature. Specifically, let S(t, z) satisfy the one-parameter group 
Equation |^ hence, PDE §§, with v(z) = %S{0,z). Then t i-> S(t,z) 
solves a separable ODE with parameter z, Equation |^. Thus, the 
canonical parametrization 7(t,x) = S(—t/2,x) itself satisfies the ODE 

d*y 

(26) — = i v ( 1 ) ) j(0,x) = x 

In particular, W(z) = iv(z)d z may be interpreted as the normal velocity 
vectorfield of the moving curve T t = {^(t,x) : x G M.}. Consequently, 
the singular foliation {T t } coincides with the family of equipotentials 
defined by the dual meromorphic differential uj = dz/iv(z): 

(27) Re[dz/iv(z)) = 

Locally, writing W = and uj = d<p, we have Re[4>] = const = C{t) 
along the equipotential curve T t . 

Proof. We begin with a more direct derivation of Equation Assume 



Equation |2l| holds. Then differentiation of the identity S(t, S(—t, z)) = 
z with respect to z and t yields two equations: S'(t, S(—t, z))S'(—t, z) = 
1, and S(t, S(— t, z)) = S'(t, S(—t, z))S(—t, z). We thus obtain the 
equation S(t, S(— t, z)) = S(—t,z)/S'(—t,z) = —2iv(z). Therefore, 
—2iv(S(t, z)) = S(t, S(—t, S(t, z))) = S(t, z) holds, hence, Equation |25|. 
In view of the previous proposition, Equation ^6| holds as well. We re- 
mark that SAz) also satisfies an ODE in z, for fixed t, obtained by 



combining Equations 22 and 25 



(28) ^ = v(S)/v(z) 

Note that iW = —v(z)d z is tangent to the foliation {T t } and Equa- 
tion ^7] follows from uj{iW) = i. (Here it is convenient to identify a 
(l,0)-complex vector X = (a+ib)d z with its real part Re[X] = ad x +bd y ; 
but one could as well denote tangent and normal vectors to {T t } by 
i?e[zVK] and i?e[W], respectively, and write the equation for equipoten- 
tials of uj = (a + i(3)dz as = i?e[u;](i2e[X]) = Re[u(X)} = aa — (3b.) 

For the last statement, we use the fact that a holomorphic differential 
uj is closed, so can be written locally as dip = <p'dz for some complex 
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potential <fi — U + iV (and the dual vector to d(f> is 8$ = \{djj —idy) = 
(l/4>')d z ). Equation |27] is the derivative of Re[4>\ = const along T t . □ 

To compare the dynamics of T t with a stationary ideal flow, we 
regard the differential uj = fdz as the starting point. For simplicity, 
we assume the local form uj = d<ft = <ft'dz, with complex potential 
(f) = U + iV defined on some planar domain. The fluid streamlines 
Im[dcf)\ = (or Irn[<p] = V = const) are interpreted as the trajectories 
of fluid particles with velocity V£7 = U x + iU y (in "scalar" notation). 
Indeed, the Cauchy-Riemann equations imply: 

a) V£7 is perpendicular to W = V x + iV y , and hence trajectories of 
V£7 are level curves of V, 

b) VU is locally sourceless and irrotational — z = VC/ is ideal, 

c) VU(z) = U x - iV x = 

Meanwhile, the equipotentials are the orthogonal curves, Re[d(f>] = 
(or Re[4>] = U = const). The latter may also be regarded as streamlines 
of the conjugate flow W associated with the conjugate differential 
given (in the holomorphic case) by *d<p = —idip. 

While the differentials, conjugates, streamlines and equipotentials 
are invariant notions, the gradient field VU implicitly depends on the 
metric induced by the coordinate z. Reverting to the operator no- 
tation for vectors, we define (Ad z ,Bd z ) = AB; then Y = <fi'(z)d z 
is ( , )-dual to d(j) = (fi'dz in the sense that, for each X = Ad z , 
(X,Y) = A(p'(z) = d(p(X). Finally, freeing ourselves from coordi- 
nates, we consider "generalized fluid flows" dual to u with respect to a 
globally defined conformal metric ( , ). In this setting, the relationship 
between the flow of T t and ideal flows may be summarized: 

Corollary 10. For a given differential uj = (ft'dz, the dual velocity 
field W = d<f) of the moving curve T t is reciprocal to the ideal flow 
Y which is ( , )-dual to uj. That is, for any ( , ) , W and Y have the 
same directions, but reciprocal magnitudes. 

Proof. First consider a local metric determined as above by coordinate 
z. Then W = = §§<9 Z = (1/0')^ = (¥U)IW{z)\ 2 )d z , while Y = 
<f)'(z)d z . Now if ( , ) is replaced by a conformally related metric, say 
((Ad z ,Bd z )) = X 2 AB, then the ((,)) -gradient, Y = X~ 2 Y satisfies 
((Y, Y)) = X~ 2 ( Y, Y) = 1/A 2 ( W, W) = W, W)). □ 

Example 11. Point vortex vs rotating line The differentials uj = dz/ z, 
*uj = dz/iz and standard planar metric yield the conjugate pair of 
ideal flows z = x + iy = (x + iy)/(x 2 + y 2 ) ('point source') and 
z = (—y + ix)j{x 2 + y 2 ) ('point vortex'); the corresponding reciprocal 
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fields x + iy and — y + ix are, respectively, velocity fields for continu- 
ous reflection through expanding circles and rotating lines. Evidently, 
such reciprocal flows need not be irrotational nor divergence-free. The 
rotating line is the special case alS/ S = —2idt of Equation |25], with so- 
lution S(t, z) = exp(—i2t + log z) = e~ 2lt z, while the expanding circle, 
S(t,z) = exp(— 2t — log z) = e~ 2t /z, solves dS/S = —2dt. By anal- 
ogy with classical terminology, the latter may be called conjugate to 
the former. (Note the expanding circle does not fit our normalization 
S(0, z) = z; the use of general base curves, in the sequel to this paper, 
will facilitate the discussion of such conjugate pairs of solutions.) 

The rotating line S(t, z) = e~ 2lt z is a well-defined global solution, 
despite the fact that separation of variables produced the multival- 
ued function log z in the intermediate expression. The presence of the 
stationary point is responsible for this slightly subtle behavior. In a 
neighborhood of a non- stationary point, a solution to Equation |25] is 
locally conjugate to the translation solution S(t, z) = z + it (the case 
v(S) = -1/2 of Equation |2J). 

Indeed, if v(x) is a real analytic function which does not vanish at 
Xq G R, antidifferentiation yields an open map F(z) = J xq _2v(s) w ^ n 
local inverse h(z) = F' 1 ^). Then for small \t\ and for z near xq, the 
solution to Equation ^5] may be expressed as: 



(29) S(t,z) = F~\it + F(z)) = hoT t oh-\z), T t {z)=z + it 



(For general initial conditions — 5(0, z) no longer the identity — conjugation 
gets replaced by Hermitian conjugation as in Equation [13].) 

On the other hand, suppose g(z ) = —2iv(z ) = 0. Then §j.S(t, Zo) = 
g(zo)S'(t, zo) = 0, so S(t, Zq) = S(0, Zq) = z§\ for all t, Zo is a fixed point 
of z i — ► S(t,z). When zq = xq is real, x = S(t, Xo), so xo lies on each 
T t — i.e., xo is a stationary point. Such xq is a pivot point if v'{xq) ^ 0, 
and a point of tangency if v'(xq) = 0, as explained in the following 

Proposition 12. Let S(t, z) satisfy Equation let xo be a real zero 
ofv(z), and consider the residue X = Res(l/v(z);x ). Then the rate 
of angular rotation of the tangent line to T t at the stationary point Xq 
is constant in time, and is given by: 



(30) 



de 

~dt 



v'(x ) 



A 



1 
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The second equality holds when x is a simple zero; otherwise, T t re- 
mains tangent to the real axis at xq, while the curvature and arclength- 
derivative of curvature of T t change at constant rates, as does the Kas- 
ner invariant, when A ^ 0: 

(31) ^ = (v"(x )f/v"'(x ) = ± 

In particular, the rates % and ^ are local conformal invariants for the 
continuous reflection equations. 

Proof. Since xq lies on the curve T t , for all t, we may define a function 
of time a{t) = S'(t, x ), which Equation || allows us to write as a{t) = 
e- 2 ^K Setting z = x in ^S{t,z) = g'(z)S'(t, z) + g(z)S"{t, z), we 
see that a satisfies the ODE: % = —2iua, which is equivalent to fl30|). 

Next, when g(x ) = g'(x ) = 0, we have S'(t,x ) = 1 since x is a 
point of tangency. Therefore, evaluating ^- t S = g" S' + 2g' S" +gS"' at 
z = xo, gives S"(t, xo) = g"(x ). The t-derivative of Equation || at z = 
xq is thus = ^g"(xo) = v"(xq). Similarly, S"'(t,xo) = g"'(xo) + 

3g"(xo)S"(t,x ), and the t-derivative of Equation |9] at z = x gives 
= i(S"'(t,x )-3S"{t,x )S"{t,x )) = ig"'(x )=v»'{x ). Since 
k(0,x o ) = k s (0,x ) = 0, we find that n(t, x ) = tv"(x ) and K s (t,x ) = 
tv'"(x ), and arrive at JC(t,x ) = t(v"(xo)) 2 /v"'(xq). The expression 
involving A is readily verified; a more comprehensive discussion of the 
significance of the residue A is given in [|Calini-Langerj1. □ 



Example 13. Continuous Reflection of Circles 

Letting u, A, and B of Example |l] now depend on t, (p0| ) becomes: 

^ ^ = g{z) = i((Auj-Au)z 2 +i(uu!-uJiv+AB-AB)z+(uB-ujB)), 

where g(z) = i{c<iz 2 + c\Z + Cq) has constant coefficients. The implied 
ODE system may be used to show that the (lifted) solutions correspond 
to intersections of the hyperboloid with two-dimensional planes 

through the origin (fitting a standard description of the geodesies in 
space forms). 

More to the point, we now illustrate the distinct types of solutions 
S(t, z) resulting from separation of variables in Equation |25|. Using 
S(0, z) = z (w(0) = 1, A(0) = B(0) = 0), we find that T Id A 3 = {g(z) = 
—2iv(z) = i(a 2 z 2 + a\Z + ao) : e R}. Consider the following three 
examples, in which v(z) has, respectively, a pair of complex zeroes, a 
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2 z cosh t + i sinh t 



-iz sinh t + cosh t 
z cos t + i sin t 



pair of real zeroes, and a double real zero: 
dS/(l + S" 2 ) = idt : S(t, z) = tan(it + arctanz) : 

dS/(l — S 2 ) = idt : S(t, z) = tanhfzt + tanh -1 z 

> y 1 y ' v ' z^sint + cost 

dS/S 2 = idt : S(t, z) = -l/(it - l/z) = z/(-itz + 1) 

We have normalized the expressions on the right-hand side to apply a 
standard criterion: if S = satisfies det(S) = ad — be = 1 and r = 
tr(S) = a + d e E, then S is hyperbolic, elliptic, or parabolic, depending 
on whether tr(S) > 2, tr(S) < 2, or tr(S) = 2. (The loxodromic non- 
hyperbolic case, tr(S') 7^ rea/ does not occur for Schwarz functions.) In 
fact, S(t, z) describes a pencil of circles associated with (respectively) 
two complex fixed points, two real fixed points, or one real fixed point. 
In the hyperbolic case there are no stationary points on T t , in the 
elliptic case there are two pivot points, and in the parabolic case there 
is a point of tangency — see Figure [I]. (The elliptic and parabolic cases 
allow Xo = 00, as in the rotating/translating line solutions already 
discussed.) 

The corresponding geodesies in Q\(l) have initial vectors g which are, 
respectively, timelike, spacelike, or lightlike, as vectors in 1^; in fact, 
the causal type of g(z) = i(a 2 z 2 + a,\Z + a ) G T Id A 3 ~ 7 1 n ) o ) o,o)Qi(l) 
is determined by the discriminant 4(g,g) = a\ — 4a2do- (We note also 
that r(t) satisfies the initial value problem: r tt = —(g,g,)r, r(0) = 
2, Ti(0) = 0, directly relating the causal type of g to the type of pencil 
given by S(t, z).) 




Figure 1. Iterated reflection of circles: the hyperbolic, 
parabolic and elliptic cases. 
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In the above example, solutions S(t, z) are globally defined, and all 
curves T t are regular. We now illustrate how continuous reflection may 
be meaningfully continued through a singular curve T to . 

Example 14. Confocal Ellipses 

Formally solving dS/ y/1 — S 2 = —idt, S(0, z) = z gives S(t, z) = 
cos(it + cos -1 z) = z cosh t — iyl — z 2 sinht. The ODE requires inter- 
pretation due to multivaluedness, which also doesn't quite resolve itself 
in the final expression for S(t, z) (as in earlier examples). Observe that 
the images of the horizontal line segments 7i t = {x + it : < x < 2tt} 
under the map h(z) = cosz = cos x cosh y — zsinrcsinhy are (once- 
covered) confocal ellipses with focii at ±1 (and the images of vertical 
lines are the confocal hyperbolas orthogonal to the ellipses) . This sug- 
gests S(t, z) = hoT t ohr l [z) (see Equation ^) should be interpreted as 
an evolution through confocal ellipses, which degenerate at t = to the 
line segment [—1, 1] — see Figure Q the line segment may be thought of 
as doubly traversed, and motion through t = is better visualized if 
the ellipses are oppositely oriented on either side of t = 0. 

To realize this description, yl — z 2 (likewise yl — S 2 ) needs to be 
defined on a Riemann surface constructed out of two copies D± of the 
slit domain D = C — [—1, 1], glued along the branch cuts [—1, 1] in the 
usual way. For t > 0, z and T t lie on the sheet D + , and for t < 0, 
z and F t lie on the sheet D_. By adding two copies of oo ("north 
and south poles"), the surface may be completed to become a copy 
of the Riemann sphere S 2 , on which T t becomes an upward- moving 
circle, coinciding with the equator at t = 0, and approaching ±oo 
asymptotically as t — ► ±oo. Meanwhile, the conjugate solution, which 
sweeps out the orthogonal family of hyperbolas, metamorphoses to a 
circle pivoting about north and south poles on S 2 . Up to conformal 
equivalence, we have merely reconstructed two of the three cases of 
the previous example — the ellipses being hyperbolic, the hyperbolas 
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elliptic! Nevertheless, this example illustrates how solutions may be 
extended past branch points via introduction of a Riemann surface. In 
the next section, we will see that there is a simpler (though perhaps 
less revealing) way to handle singularities of this type. 



4. The affine connection on A and geodesics 



Writing Equation |20| as S/S' = g(z) and differentiating yields: 



(32) S-p' = 

The main goal of the present section is to interpret this PDE as the 
geodesic equation on A, the space of analytic curves. Before adding 
geometric structure to A, however, we make a direct argument casting 
Equation 1313 as the continuous limit of Davis iteration. 



The following argument (which served as the starting point for our 
investigation), may be regarded as the "physical derivation" of the 
second order PDE for continuous reflection. To start, we note that 
Davis iteration may be written 



(33) 5**0+2* = S to+t o Sto o S, 



to+t 



where to and t are taken to be integers (in fact, this equation follows 
from Equation |T7| by elementary symmetric space identities given in 
the appendix). 

However, by scaling time, to and t ought to be allowed to assume 
any real values; indeed, from now on, we simply regard to an d t as con- 
tinuous variables. On the other hand, to derive our second order PDE, 



we will assume only that Equation [33| holds to second order in t, for 
each t . Thus, we differentiate twice with respect t at t — 0, and apply 
first and second order identities satisfied by Schwarz functions (such 
identities will be useful later on as well). It is notationally convenient 
to set to = for computation, but it will be clear that we may reinsert 
to in the resulting equation. First we take of the right-hand side: 

^S(t, S(0, S(t, z))) = S(t, S(0, S(t, z)))+S'(t, 5(0, S(t, z)))S'(0, S(t, z))S(t, z) 
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Taking second derivative leads to six terms, two of which combine: 

dt 2 ' 



^RHS = S(t,S(Q,S(t,z))) 



+ 2S'(t, 5(0, S{t, z)))S'(0, S(t, z))S{t, z) 

+ 5"(t, 5(0, S(t, z)))[S'(0, S(t, z))S(t, z)f 

+ S'(t, 5(0, S(t, z)))S"(0, S(t, z))[S(t, z)f 

+ S'(t, 5(0, S(t, z)))S'{0, S(t, z))S(t, z). 

Note that the identity S(t, S(t, z)) = z implies S'(t, S(t, z)) = 1/S'(t, z) 
and also S"(t, S(t, z)) = -S"(t, z)/[S'(t, z)] 3 . Applying the three iden- 
tities at t = gives 

—RHS\ t=0 = 5(0, z) + 25'(0, z)S(0, z)/S'{0, z) + 5"(0, z)[S{0, z)/S% z)f 
ot z 

+ 5(0, z)(-5"(0, ^)5 2 (0, z)/5 /3 (0, z)) + 5(0, z)S(0, z)/S'{0, z) 
= 2S(0,z) + 2S'(0,z)S(0,z)/S'(0,z) 

Equating this last expression to -^LHS\t=o = 45(0,2;) and rearrang- 
ing, we obtain Equation 32 at to — 0; however, the derivation is ac- 
tually valid for any time to- Thus, we have recovered Equation [32|, 
and by antidifferentiation, Equation |2(| What we have gained by the 
second order interpretation of continuous reflection is comparability to 
some familiar classes of physical and geometrical systems; in particular, 
continuous reflection may be regarded as a local, autonomous process 
(which does not require us to unnaturally postulate exact interpolation 
of a discrete process). 

In the remainder of this section we consider time-dependent curves 
w(t,x) obtained from non-singular holomorphic functions z i— > w(t,z) 
defined for z in a neighborhood of the real axis, and we consider also 
the corresponding time-dependent Schwarz functions defined by the 
identity w(t,x) = S(t,w(t,x)). In fact, we assume w to be invertible, 
as in the previous section, and freely apply the substitution z = w(t,x), 
x = u> -1 (t, z) and analytic continuation as necessary, e.g., to obtain the 
representation of Schwarz functions: S(t,z) = w(t,w~ l (t,z)). Here we 
abuse notation by writing w~ l (t,z) for the inverse of the map z i— ► 
w(t,z). (For the rest of this section, we reserve subscripts for partial 
differentiation, so the notation w t ~ 1 (z) is unavailable. Note also the 
letter l w' has replaced '7', which by now is reserved for the canonical 
parametrization along geodesies.) 

A (formal) tangent vector to A = {T} ~ {5} at a Schwarz function 
S(z) = 5(0, z) is understood to mean a complex analytic function along 
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{z : z = S(z)} coming from an infinitesimal variation of Schwarz 

at 1 



functions, X = §fS(t, z)\ t=0 . Given tangent vectors X,Y at S, we will 



need to consider quantities such as X z , X 2 , X Z Y, defined by the usual 
operations of differentiation and pointwise multiplication of functions. 
Taking further advantage of the inclusion of A in a linear space TC 
of holomorphic functions (defined on a suitable domain U D T), we 
form also VxY, the covariant derivative of Y with respect to X in the 
ambient space 7i. Here, X and Y are vectorfields on A, extended to a 
neighborhood in H to define VxY. If S(r, t, z) is a two-parameter vari- 
ation of Schwarz functions, and X, Y are extensions of X = -§pS(r, t, z) 

and Y = §^S(r, t, z), then we have VxY = -,^S(r, t, z). 

The special combinations of the above operations on X, Y yielding 
again tangent vectorfields to A may be identified using the following 

Lemma 15. Let V be an analytic curve with Schwarz function S(z). 
Then the tangent space to A at S is given by: 



T S A = {X G H : X(z) = a(z) - S z (z)a(z) along T, for some a G H} 
= {X G H : X 2 /S z < along T} 

Proof. Care must be taken here to distinguish between conjugates of 
numbers and functions! Consider S(t,z) with 5(0,2:) = S(z). Differ- 
entiating w(t,x) = S(t,w(t,x)), we obtain S t (t,w(t,x)) = w t (t,x) — 
S z (t,w(t,x))w t (t,x). Letting z = w(x) = w(0,x), we obtain X(z) = 
S t {0,z) = w t (0, nr^(z))-S z (z)wtiP , w^jz)). Nowseta(^) = w t {0, w~ l (z)), 
and note that a(z) = w t {Q-,w~ l (z)) = w t (0, w~ 1 (z)), for z G T; thus, 
X(z) = a(z) — S z (z)at(z) along T. Observe that a(z) is an arbitrary 
analytic function along T, so the first description of T5A follows. 

Using the same representation of S and X in terms of w, and using 
the fact that S z (t,w(t,x)) = w x (t,x)/w x (t,x), we have 

X 2 , ( \ \ (w x w t - w x w t ) 2 f . 2 \2 . n 

— (w{x)) = = U= = {-2iRe[iw x w t \/\w x \) = -4{n,w t ) <0 

O z VJ X U) X 

For the normal variation, (h,w t ) = Re[iw x w t ]/\w x \, we use the posi- 
tively oriented orthonormal frame t = w x /\w x \, h — it along w(t,x). 
Note that the normal variation corresponds to an arbitrary analytic 
function along T. □ 

Theorem 3. a) Equation [51] for continuous Schwarzian reflection is 
the geodesic equation Vs t S t = for a (formal) symmetric affine con- 
nection on A, defined by: 

(34) V X Y = V X Y - ^{X Z Y + Y Z X) 

Zo z 
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b) In fact, V is the canonical connection on the (format) symmetric 
space A, determined by the multiplication dl2|); accordingly, V is con- 
formally invariant, with covariant- constant Riemann tensor. 

Proof. Using Vs t S t = -j^S(t,z), Equation ^ follows at once from 
Vs t Sf = 0. The main content of a), therefore, is the interpretation 
of V as a connection. The latter is included in b) which will be treated 
in Appendix II; however, it is useful to include a direct computational 
argument for a) here. 

Given our use of the embedding A C H, the main point we need to 
verify is that VjcF is in fact tangent to A. We introduce a notational 
shorthand suggested by the lemma: for a given curve w(x), an analytic 
function g(x) can be turned into a tangent vector X(z) = g(w~ 1 (z)) to 
A at S(z) = w(w^ 1 (z)), by the operator g(x) = g(x) — S z (w(x))g(x). 
(The formula X(z) = g(w~ 1 (z)) is valid only along T, though X(z) 
does indeed represent an element of TC, by analytic continuation). 

Now consider a two-parameter family of non-singular analytic curves 
w = w(r, t, x) and the corresponding two-parameter variation of Schwarz 
functions satisfying S(r,t,w(r,t,x)) = w(r,t,x). We compute suc- 
cessive derivatives of the latter, suppressing all arguments: w — S, 
w x = S z w x , w xx = S zz w1 + S z w xx , w r — S r + S z w r , w rx = 
S rz W x + S zz W x W r + S z W rx , W rt = S r t + S rz w t + S zt w r + S zz w t w r + S z w rt 
(and similar formulas for Wt and Wt x ). Rearranging, and applying the 
vector operator to g(x) = w xx (r, t, x), g(x) = w rx (r, t, x), etc., one gets: 

w„S zz = w xx 



x 



w x S rz — w x w rx w r w xx 
w 2 x S tz = w x w tx - w t w xx 

wlSrt = wlw rt - W t W x W rx - W r W x W tx + W r W t W, 



XX 



These formulas are substituted into 2w x w1V s r St = 2vJ x wlS rt —w x (S r Stz+ 
StS rz ) (obtained by setting X = S r , Y = S t in Equation [34|) , and coeffi- 
cients of the second order terms w xx , w rx , w tx , w rt are collected. After a 
key cancellation and subsequent division by w x , these coefficients turn 
out to be real: 2w x w x 'V s r S t = 2w x w x w rt — (w x w t + w x w t )w rx — (w x w r + 
w x w r )w tx + {w r w t + w r w t )w xx . For future reference, we re-express this: 

(35) V Sr S t = w rt + A^(K, w t )w xx - (w x , w t )w rx - (w x , w r )w tx ) 

\Wx\ 

In particular, after precomposing with w" 1 (z), the right-hand-side 
has the form of a tangent vector to A at S. Finally, replacing S r , S t 
with general vectorfields X, Y on A, one obtains (X, Y) i— > Vx7, 
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a bilinear operation satisfying the required further properties for an 
affine connection: V fX Y = fV X Y and V x fY = fV X Y + X(f)Y, for 
/ : A -> R. □ 

Corollary 16. S(t,z) = w(t, w _1 (t, z)) represents a geodesic Y t in A 
if and only if the parametrization w(t,x) of Ft satisfies Im[Q[w]] = 0, 
where Q is the quartic second order partial differential operator 

(36) Q[w] = w x w 2 x w tt + (w x w x w t + w 2 x w t )w xt + w x w t w t w xx 

In the case of normal motion, (w x ,w t ) = 0, the equation Im[Q[w]] = 
may be replaced by the system 

(37) Im[w x (\w x \ 2 w u + \w t \ 2 w xx )} = 0, w t w x + w t w x = 
Proof. Setting r = t in Equation [3^, we obtain the geodesic equation: 

(38) = V St S t = w a + T—lo(( w ti w t)Wxx - 2(w x ,w t )w tx ) 

\Wx\ 2 

One easily verifies that w 2 w x V ' s t St = —2Im[Q[w]], and the corollary 
follows. To relate the corollary more directly to Equation [3^, we note 
that w x Vs t St = w 2 (S x S t t — S t S xt ). That the geodesic equation may 
be written in the form Im[* * *] = is not surprising, given that T t 
is an evolving unparametrized curve, so the above equation is missing 
"half" of the information required to govern w(x,t). Up to initial 
parametrization, this information is provided by the normal motion 
requirement in the second statement of the corollary. □ 

Remark 17. It is evident that the equation Jm[Q[w]] = must be in- 
variant under conformal transformation and reparametrization. Specif- 
ically, one may directly verify that if w v (t,x) = <p(w(t,x)), then 

Im[Q[w v }} = \\^\\ 4 Im{Q[w}}. 

Also, if o~(t, x) is a (time-dependent) diffeomorphism of the parameter 
domain /cl, and w a (t, x) = w(t, a(t, x)), then 

Im[Q[w a }} = allm[Q[w}}. 

We remark that Q[w] itself does not satisfy either invariance property; 
in particular, i?e[Q[/uv]] turns out to be a rather lengthy expression 
(for which we have not found any interpretation). 



Corollary 18. Let the family of Schwarz functions S(t, z) satisfy the 
one-parameter group Equation f|7], for z near the real axis. Then the 
canonical parametrization j(t,x) = S(—t/2,x) satisfies Equation 
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Corollary 19. Any holomorphic function f(z) may be interpreted as a 
geodesic in A, in the sense that w(t, x) = f(x + ict) solves Equation 
for c a real constant. 

We emphasize that the converse is far from true; namely, w(t, x) = 
f(x + ict) cannot represent a solution to the geodesic equation in 
a domain including any stationary points — which, as we have seen, 
are important elements of geodesic behavior. On the other hand, 
the corollary extends our original notion of continuous reflection by 
slight of hand: the equation = im[Q[w]] is perfectly meaningful and 
well-behaved where w x vanishes — where the actual geodesic equation 
= Vs t St = —2Im[Q[w]]/wlw x is not a priori meaningful. 

Now let us consider what qualitative behavior may be exhibited by 
such singular geodesies. For convenience of working with closed curves, 
we will write instead w(t,9) = h(e l9+t ); in other words, we transform 
the expanding circle solution Ct = {e t0+t G C : < 9 < 2n} by a 
holomorphic (or meromorphic) function h. For example, with h(z) = 
^(z + -), we get w(t,9) = cos(6> — it) = cos^cosht + isin^sinht— 
precisely the family of confocal ellipses described in Example [14] (with 
orientation flipping at the singular curve w(0,9) = cos 9). 

In general, consider the discrete set of exceptional t- values {tj} at 
which Ct encounters a singularity of h (for the above ellipses, just 
to = 0). Between consecutive exceptional values, w(t,9) parametrizes 
a regular homotopy of closed curves Y t = h(C t ); that is, for tj < 
t < tj + \, < 9 < 2tc, w(t,9) has a non-vanishing and continu- 
ously varying tangent w'(t, 9). In particular, the rotation index I\Tt]~ 
the degree of the unit tangent to T t as a map of between circles, 
e ie+t ^ #)| — is constant on (tj,tj+x). At each tj, on the 

other hand, the index I[T t ] jumps: e.g., a simple critical point (where 
h' '(e j6, ° + ^) = 0, h"(e ie ° +tj ) 7^ 0) creates a new loop just after T t develops 
a cusp at h(e t9o+t j) — see Figure^; a pole of h' destroys a loop (or loops) 
by a breaking/reconnecting transition — see the Celtic Cross in Figure |j. 
Such conclusions are based on 

Proposition 20. Let h(z) be meromorphic on the open disc D = {z 6 
C : \z\ < r}, and analytic and nonsingular on the circle C = 3D. Then 
the rotation index of the image curve h(C) is given by 

(39) I[h{C)\ = Z[h'\ - P[h'} + 1 

where Z[f] and P[f] denote, respectively, the number of zeros and poles 
of a meromorphic function f(z) in D, counting multiplicity. 

Proof. The curve h(C) parametrized by w(9) = h(re ld ), < 9 < 2tt, 
has tangent vector w'{9) = ire 10 h' \re ld ) . Setting g(z) = izh'(z), it 
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follows that the degree of the unit tangent to h(C) is the winding 
number of g(C) about the origin. Applying the argument principle to 
the curve C and the function g(z), we find that 

I[h(C)] = ^A G arg[g(z)\ = Z[g] - P[g] = Z[h'\ + 1 - P[h'\ 

We remark that the curve C may be replaced by an arbitrary simple 
closed curve, as in the argument principle, but such a generalization is 
of no importance to us presently. □ 



15 20 25 



Figure 3. Reflection through cusps (h(z) = z 2 -Az + A } 
and h(z) =A/(z 4 + Az)) 



Example 21. Singular geodesies via rational functions 

For a rational function h(z) = p(z)/q(z), viewed as a map into the 
Riemann sphere, the breaking/reconnecting transitions are replaced by 
regular homotopies in which the curve w (t, 8) = h(e ie+t ) moves past oo 
(= north pole). Note that the poles of h' = (qp' — pq')/q 2 are of even 
order; at a pole of order 2, e.g., a positive loop is converted to a negative 
loop, resulting in a net change AI = — 2 (for the stereographic image 
in the plane). Here one should recall the Whitney- Graustein Theorem 
and its analogue for curves on S 2 : while regular homotopy types of 
planar curves are classified by J, regular curves on S 2 are (effectively) 
classified by Imod 2 (see | [5 maid] , §7). To summarize: w(t,6) = h(e t9+t ) 
parametrizes a regular homotopy of curves T t C S 2 , deleting at most 
deg(p) + deg(q) — 1 values of t, and the regular homotopy type of 
T t changes where (qp' — pq') has a zero of odd order. With these 
considerations in mind, one may construct singular geodesies having 
more or less prescribed cusping and looping behavior. A few examples 
are pictured in Figures H, |. 
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5. Appendix I: symmetric space formalism 

In this section we briefly review enough of the basic definitions and 
algebraic constructions of symmetric spaces — as developed in Symmet- 
ric Spaces I, by Ottmar Loos — to explain how the formal structure of 
Schwarzian reflection fits into the theory. 

Definition 22. (Loos): A symmetric space is a manifold Ai with a 
differentiable (non- associative) multiplication fi(p,q) — p- q satisfying: 

1. p ■ p = p 

2. p - (p • q) = q 

3. p ■ (q ■ r) — {p • q) • {p ■ r) 

4. Every p has a neighborhood U such that p ■ q = q implies p = q 
for all q G U . 

(The definition of a locally symmetric space is similar, with /i replaced 
by a germ of a multiplication about each point p G Ai.) 

To interpret (l)-(4), consider the mapping q \— > s p (q) = p(p, q), called 
left multiplication by p, or the geodesic symmetry of Ai with respect to 
p. Then (2) and (3) may be written s^ = Id and s p (q-r) = (s p q) ■ (s p r). 
In other words, s p is an involutive automorphism of Ai; further (1) and 
(4) state that p is an isolated fixed point of s p . 

More standardly, one views the multiplication /i on a symmetric 
space M. as a byproduct of existing geometric structure (say, Rie- 
mannian or affine), in which case s p is an isometry defined by re- 
versing geodesies through p — hence the name. In any case, one pro- 
ceeds to generate the group of displacements, G(M), from the auto- 
morphisms/isometries d PA = s p s g , using all pairs of points p, q G M.. 
This normal subgroup of the group of symmetric space automorphisms, 
Aut(M.), should be thought of as "large"; e.g., G(A4) is transitive on 
connected Ai, and G(A^) = Aut(Ai) , the identity component of the 
automorphism group, for semisimple Ai. 

In the "algebraic" approach to symmetric spaces, it takes some work 
to develop all the familiar geometric structure out of the multiplication 
alone. On the other hand, one arrives very quickly at the large group of 
symmetries G(Ai). Thus, the approach is not far removed from Klein's 
idea, to define a geometry by its symmetries. 

We now list several (closely related) symmetric space constructions 
(and refer the reader to |[Loos|| , for more details and examples). G will 
always denote a connected Lie group. Generally, axioms (l)-(3) are 
easily verified by direct computation, while (4) usually requires a bit 
more work. 
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Lie groups with symmetric space product (G,n): For any G, 
one may introduce the symmetric space multiplication 

(40) fi(g, h) = g ■ h = ghT l g 

For example, if V is any vector space, regarded as an additive group 
G = V, the product x ■ y = s x y = 2x — y agrees with the familiar 
reflection ofy in x. The spaces (G, /x) play an important role as ambient 
spaces into which other symmetric spaces may be canonically embedded 
as symmetric subspaces (as explained below). 

We note that for finite dimensional groups, (4) follows easily from the 
existence and basic properties of the exponential map. Some infinite 
dimensional groups — in particular, loop groups — pose no difficulty in 
this regard. On the other hand, for diffeomorphism groups, the bad 
behavior of the exponential map is problematic, and leads to interesting 
subtleties. 

Symmetric subsets G a : Assume there exists a : G — > G an 
involutive automorphism: o~{gh) = o~(g)o~(h), and cr 2 (g) = g. Then the 
subset G a C G of symmetric elements (or transvections) is defined by: 

(41) G a = {a(g)g- 1 : g G G} 

In other words, G a = C{G) is the image of the Cartan submersion: 

(42) C : G - G a , C(g) = a(g)g- 1 

Note that an element p G G a satisfies p^ 1 = o~(p); it can be shown that 
G a is in fact the identity component of {p G G : p^ 1 = o~(p)}. 
Now G a is closed under the multiplication 

(43) p • q — pq~ x p = pcr(q)p, 

which agrees with the one given above on G itself. Thus, G a is a 
symmetric subspace of G. Further, there is a natural left action of G 
on G a , given by 

(44) X(g)(p) = a(g)pg' 1 

In fact, A defines a homomorphism of G into the group of symmetric 
space automorphisms Aut(G a ). 

Symmetric homogeneous spaces G/K: Assume again an in- 
volutive automorphism o~(g) on G. Let G a be the fixed point set 
G a = {g G G : a(g) = g}, and G a its identity component. For K 
a (necessarily) closed subgroup of G satisfying G a C K C G a , we con- 
sider the quotient G/K, with elements [g] = gK. Using the fact that 
a is an automorphism, one checks that the formula 

(45) \g\ ■ [h] = [ga{g)-'a{h)\ 
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gives a well-defined symmetric space multiplication on G/K. 

Conversely, the homogeneous space structure of an abstract sym- 
metric space may be recovered from the multiplication, via the group 
of displacements. Specifically, assuming a base point o in M. has 
been fixed, let H C G{M) be the isotropy subgroup of o. Then 
a(g) = SogS' 1 = S gS defines an involution of G(A4) (recall, G(Ai) is 
normal in Aut(M)), G a Q C H C G° ', and M is isomorphic to G(M)/H, 
with the symmetric space structure on the quotient defined as above 
(see ||Loosfl , Theorem 3.1). 

Next, we describe the relationship to the previous symmetric space 
construction. We use the same letter C to denote the Cartan immer- 
sion: 

(46) C : G/K - G a , C{[g\) = a{g)g- 1 

Then, as is easily verified, C gives a well-defined symmetric space ho- 
mo morphism C : G/K — > G a , which becomes an isomorphism in the 
special case K = G G of the symmetric homogeneous space G/K. Cor- 
responding to A, the left action of G on G/K is given by left translation, 
defining a homomorphism of G into the group of symmetric space au- 
tomorphisms Aut(G/K): 

(47) r(g)([h}) = [gh] 

More explicitly, C intertwines the two actions: Cr(g)([h}) = X(g)C([h\). 

We note that for M. = G/K, the group of displacements, G(A4), 
is precisely the group generated by r{G a ) = r(C(G/K)). In fact, 
using the natural base point, o = [e] = K G G/K, consider the 
quadratic representation, Q : M. — > G(M), defined by Q{p) = s p s Q . 
Then we have Q([g])([h}) = s [g] s [h] = [g] ■ ([e] • [h]) = [g] ■ [a(h)\ = 
[ga(g)- l h] = r(C([a(g)]))([h}), so Q([g]) = r(C([a(g)})). On the other 
hand, Q{p)Q{q) 1 = s p s (s q s ) 1 = s p s q 1 = s p s q . Thus, G{M) is 
generated by r{G a ) = r{C{G/K)). 

Symmetric Quadric Surfaces Q™(pi): Let be the n- dimensional 
pseudo-Euclidean space with metric 

(48) (x, y) = x 1 y 1 H h x n _ s y n _ s - x n _ s+1 y n _ s+1 x n y n 

Then the quadric surface Q™(a) = {x G : (x, x) = a} is closed 

under the product: 

(49) x . y = 2 ^l4x-y 

In this one example, (3) requires a bit more work, whereas (4) follows 
directly from the fact that x ■ y = y implies y = ±x. 
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Among the standard spaces obtained this way are the (Riemann- 
ian) spheres and hyperbolic spaces, Q\(— 1), as well as Lorentzian and 
pseudo-Riemannian space forms. Of course, the spaces Q™ are related 
to the classical groups, hence also to the other constructions listed 
above. On the other hand, there are advantages to the present direct 
approach. For instance, note how this example generalizes without dif- 
ficulty to infinite dimensions. Namely, let (Hi, ( , )j) be Hilbert spaces, 
for i = 1,2. Consider the vector space H = Hi © H2, with indefinite 
metric (x, y) = (xi, 2/1)1 — (X2, 2/2)2- Then the above multiplication on 
the quadric hypersurface QHi,H 2 ( a ) = {x G 7^ : (x, x) = a} satisfies 
axioms (l)-(4), by the identical arguments. 

Discrete symmetric spaces: Specializing the above definition to 
O-dimensional manifolds, one has simply a set M. with multiplication 
fi satisfying (l)-(3). For instance, we will regard the integers, Zcl, 
as a symmetric space with multiplication m ■ n = 2m — n. For a more 
interesting example, take a finite set of circles on the Riemann sphere, 
and iteratively generate all possible reflections (products) between pairs 
of circles. (Depending on the choice of circles, the resulting discrete 
symmetric space may or may not be discrete as a subspace of all circles). 

In a "large subspace" limit, one might regard such spaces as discrete 
approximations to continuous — even infinite dimensional — symmetric 
spaces; alternatively, the full space may be regarded as a formal sym- 
metric space, by assigning it the discrete topology. 

The following close relative of the key example A 3 C A illustrates 
several of the constructions mentioned above. 

Example 23. Oriented circles in the Riemann sphere 

On the group of Mobius transformations of the Riemann sphere, 
G ~ PSL(2, C) ~ SL(2, C)/{±1}, we consider the involution a(g) = g, 
well-defined by complex conjugation of matrices in SX(2,C), and the 
subgroup K = G u a ~ PSL(2, M) of Mobius transformations preserving 
the real axis and its orientation. We obtain the symmetric space of 
oriented circles in the Riemann sphere, M 3 = PSL(2, C)/PSL(2, R). 
A given Mobius transformation g parametrizes a circle with real vari- 
able x, and [g] = gK represents an equivalence class consisting of 
reparametrizations of g by elements of PSL(2,M). 

We note also that M 3 is naturally identified with the Lorentzian man- 
ifold Qx(l) (see |pecil|| , p. 15), thus describing the metric structure on 
M 3 . Further, the product (|49|) on Q 3 (l) induces a well-defined symmet- 
ric space multiplication on the quotient A 3 = M 3 /{±} ~ {unoriented 
circles}; in fact, the product of circle p = {±x} with circle q = {±y} 
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is none other than the reflection of q in p. Note also the implications 

x • y = y y = ±x 

x • y = — y •<=>- (x, y) = •<=>- {circles x and y meet orthogonally} 

Thus, we have simple interpretations of the fourth symmetric space 
axiom for the spaces M 3 and A 3 . Incidentally, the symmetric space 
axioms appear to preclude extending the space M 3 to include point 
circles, as in Lie sphere geometry. For continuity would require x • y = 
x, for any point circle x and disjoint circle y; hence, x ■ (x ■ y) = x • x, 
so the first two axioms cannot both hold. 

Finally, we note that for M 3 , the displacements s p s q generate the 
whole Mobius group: G(M 3 ) = G = PSX(2, C) — a classical theorem. 
One may ask if a corresponding result holds for A, at least locally: 
is any holomorphic function close to the identity locally a product of 
(pairs of) Schwarzian reflections in analytic curves close to R? 

We conclude this section with a construction giving discrete ana- 
logues of geodesies in an abstract symmetric space M. — namely, the 
powers of an element p e M. relative to a base point o e M.: 

(50) p° = o, p 1 = p, p n+1 = p n ■ p n - 1 

Note that by axiom (2), the above equation may be rewritten p n ~ l = 
p n -p n+l , so that indeed negative as well as positive integer powers of p 
are inductively determined. We note also that in the case of a Lie group 
M. = G, the powers just defined agree with the usual ones, provided 
we take as base point the identity, o = e e G. 

Proposition 24. Consider the integers 7h as an additive group, with 
corresponding symmetric space structure, m-n = 2m— n. Then the map 
m i— > 7(771) = p m defines a homomorphism 7 : Z — > M. of symmetric 
spaces. In particular, for all integers m, n, k: 

/^J p^^ ' p^ — pTfi-n — p*2)7i n ^ g ptn-{-k — ' 

a) p n+2 = p k+1 ■ (p k ■ p n ) = Q(p)p n , 

Hi) 7 intertwines respective "k th -power' ', "time-reversal" and "trans- 
lation" maps in Z and M: {p n ) k = p nW = p kn ; (o ■ p) n = p°' n = p~ n ; 
the powers of q 1 = p l+k with base point q° = p k are q n = p n+k . 

Proof. We prove the second of the two equivalent statements in i) by 
induction. For a given positive integer K, assume the equation holds 
for all 771, whenever < k < K (the case K = 1 being known). 
Using axiom (3) and (both versions of) the definition of p n , we ob- 
tain p m +( K + 1 ) — p m + K . prn+K-1 _ ^pm . pm-K^ . Q^m . pm~(K-l)^ _ 
p m . (yn-K . p (m-K)+lj = p m . p (m-K)-l) = p m . p m-(K+l) ^ g0 the regult 
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holds for K + 1. One argues similarly for K < 0. Thus, i) follows by 
induction, and 7 is a homomorphism. 

The first equality in ii) now follows easily from the homomorphism 
property, and k = gives p n+2 = p 1 ■ (p° -p n ) = s p s p n = Q(p)p n (which 
is the definition of powers given in ||Loos]| ). To prove iii), set q = p k , and 
q(n) = p kn . Note q(0) = p°, q(l) = q, and q(n)-q(n—l) = p kn .p k \ n ~ l > — 
pk(n+i) _ q(^ n _|_ Thus, we make the identification q{n) = q n , i.e., 
pkn _ (p fc )"- ; the first equation in iii) — note that in Z, the k th power 
of n is rv- k > = nk. Setting k = — 1 gives the second equation. The 
third equation follows by a similar argument using, instead, q(0) = p k , 
q = q(l) = p k+1 , and q(n) = p n+k . □ 

6. Appendix II: the canonical connection on A 

We have thus far invoked formal properties of symmetric spaces for 
heuristics, while the validity of our differential equations did not depend 
on such interpretation. We now start to consider geometric construc- 
tions which are difficult to assign precise meaning to, in the absence of 
a smooth structure on (A, fi) — yet to be defined. Thus, in the context 
of unparametrized analytic curves, Schwarz functions, etc., the defi- 
nitions and computations to follow are to be understood at a purely 
formal level. 

We begin by recalling the definition of the canonical connection on 
a symmetric space (A4 , /i) : 

(51) V X Y = X(Y) + \x-Y 

Here, X and Y are vectorfields on M. (first order operators), and we 
exploit the fact that the space of differential operators V : J-'(J^i) — > 
^(.M) may be made into an algebra in two different ways. On the 
right-hand-side, the first product X(Y) is the second order operator 
X(Y)[F) = X[Y[F}} (the notation X(Y) in place of the usual XY is 
to avoid eventual confusion with complex multiplication). The second 
term depends on the product n{p,q) = p ■ q on Ai, which induces 
the JF(7V1)-Iinear (non- associative) product on the tangent algebra of 
differential operators on M. — see [Loos] for a thorough explanation. 

Here we need the bilinear form F(X,Y) = |X • Y on TAi, which 
has values in T 2 Ai, the space of second order operators on J-{M). 
Specifically, X ■ Y[F](p) = X ®Y[F o (j](p) is computed by allowing 
the vector Y to operate on q ^ F o /i(p, g), then letting X operate on 
the resulting function of p (or vice versa), and subsequently evaluating 
at q — p. It will suffice to consider pairs of vectorfields of the form 
X = g 2 , Y = coming from two-parameter variations p{r,t) with 
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p = p(0, 0), and set 

(52) X ■ Y[F\(p) = —-F(p(r,0)-p(0,t))\r=t=o 

For the particular case of interest here, the symmetric space (A4, fi), 



defined via Equations [II], f|3, comes from a group of diffeomorphisms 
and all computations are based on the chain rule. Further, the diffeo- 
morphisms are naturally included in a linear space of functions/mappings, 
and second derivatives, etc, are simpler to discuss. 

In this setting we complete the proof of Theorem [], by carrying out 
the above computations for the concrete case M. = A. First, we con- 
sider a two-parameter variation of Schwarz functions, S = S(r,t,z), 
along with one-parameter variations P = P(r, z) = S(r, 0, z) and 
Q = Q(t,z) = S(0,t,z). Using the dot notation for the symmetric 
space product, P ■ Q{r, t, z) = P(r, Q(t, P(r, z))), and denoting partial 
derivatives of S by subscipts, we obtain the following formulas (the 
first two correspond to Lemma 2.1, p. 76 in [Loos]): 

—P ■ Q(0, 0, z) = 2S r (0, 0, z), -P- Q(0, 0, z) = -S t (0, 0, z), 



<)2 -P ■ Q(0, 0, z) = —^(S rz S t + S tz S r 



drdt ^ v ' ' ' S z 



\(0,0,z) 



The derivation is similar to that of Equation [32|, and uses the chain 
rule identities: S o S = Id, S z o S = 1/S Z , S t o S = —S t /S z , S zz o S = 
-S zz /S z , and S zt o S = (S zz S t - S z S zt )/S z . 

Regarding S r (0, 0, z), and S t (0, 0, z) as elements of Ts( 0) o,2:)A, we now 



compute the second order operator S r - S t , using Equation the above 
formulas, and the embedding A C 7i. Given F : A — >• M., we extend F 
to a neighborhood in 7i, and compute 

S r ■ S t [F] = ^ t F(P.Q)\ r=t=0 = ^DF(P.Q)(^P.Q)\ r=t=0 

= [D 2 F{P ■ Q)(—P .Q,-P.Q) + DF{P ■ Q){^P ■ Q)]r=t=o 



2D 2 F(S)(S r , S t ) - DF(S)((S rz S t + S tz S r )/S z )} r -_ 



t=o 



On the other hand, letting X = S r (r,t,z), Y = S t (r,t,z) operate 
as vectorfields in the usual way, we have SVfSif-F]] = J^(J^F(S')) = 
D 2 F(S)(S r , S t ) + DF(S)(S rt ). Putting these results into the definition 
of the canonical connection, Equation |5l], we find that the second de- 
rivative terms D 2 F cancel as expected, and we are left with: 

V Sr S t [F] = DF(S rt - (S rz S t + S tz S r )/2S z ) 
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Noting that this result does not depend on the extension of F, and that 
S r ,St may be replaced by vectorfields X, Y, one recovers the earlier 



definition of the affine connection, Equation 34 



Conformal invariance may be regarded as an "automatic" conse- 
quence of the fact that X(g)S = gSg^ 1 defines a symmetric space au- 
tomorphism and hence is affine: \{g)*V ' X Y = V \( g ) tX \(g)*Y . More 
concretely, the induced action A»:Gx TM. — > TM. may be described 
as follows. By an abuse of the above notation we have, e.g., \(g)*X = 
■jj-\(g)S = -j^g o S o g^ 1 = g z o SX\ g -i. Computing second order par- 
tial derivatives of \(g)S(r, t, z) and substituting into Equation [34|, one 
may verify: Vx( g )*xHg)*Y = g z o S(S rt - (S rz S t + S tz S r ) /2S z )\ g -i = 
g z o SV x Y\ g -i = \(g)*V x Y. 

The Riemann tensor, R(X,Y)Z = V X V Y Z - V y V x Z - V [X , Y ]Z, 
of the canonical connection on a symmetric space is always covariant- 
constant: = (V W R)(X,Y,Z) = V W {R {X,Y) Z) - R{V W X,Y)Z - 
R(X, V W Y)Z - R(X,Y)V W Z (see, e.g., flDoosfl , Corollary 1, p.84, for 
the abstract symmetric space result). Here we sketch the argument 
for the present context. We use the fact that the geodesic symmetries 
s p are also symmetric space automorphisms, hence affine: in terms 
of the induced maps s p *X — p ■ X, the bilinear form T(X, Y) = 
|X • Y satisfies p ■ T(X,Y) = T(p ■ X,p ■ Y). Now, the action of 
the induced map on tangent vectors was already given above (mod- 
ulo notation): p ■ Y — —Y. The upshot is that one is able to write: 
- (V W R) (X, Y,Z) —p- (V W R) (X, Y, Z) = {V P . W R) (p ■ X,p-Y,p- Z) = 
{V- W R){-X,-Y,-Z) = (V W R)(X,Y,Z). Since this holds for all 
W, X, Y, Z, we have VR = 0. 
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Figure 4. Gallery of singular geodesies without station- 
ary points: Weave, Celtic Cross (h(z) = 5/(z 5 — 5z)); 
Turtle (h(z) = 6(z - 3)(z + 2 - i)(z + 3i)/(z 6 -6*)); 
Bicycle Race (h(z) = 6(z-2)(z+ 1 -i)(z + i)/(z 6 -6z)) 



